TIT/HEP-581 
TU-813 
UT-08-06 
May, 2008 



Intersecting Solitons 
Amoeba and Tropical Geometry 

Toshiaki Fujimori^, Muneto Nitta^, 
J2 ' Kazutoshi Ohta'^, Norisuke Sakai"^ and Masahito Yamazaki^ 

00 ^ 

Department of Physics, Tokyo Institute of Technology 
Tokyo 152-8551, JAPAN 

p i' '^Department of Physics, Keio University, Hiyoshi, Yokohama, Kanagawa 223-8521, JAPAN 
^ '. 

^ , ^Department of Physics, Tohoku University, Sendai 980-8578, JAPAN 



OO 
O 
O 
(N 



^Department of Mathematics, Tokyo Woman's Christian University, Tokyo 167-8585, JAPAN 
> ■ 

T^-j- ■ Department of Physics, University of Tokyo, Tokyo 113-0033, JAPAN 

in 

O ' Abstract 

OO ■ 

O ' We study generic intersection (or web) of vortices with instantons inside, which is a 1/4 

BPS state in the Higgs phase of five-dimensional Af = 1 supersymmetric U{Nc) gauge theory on 
^ ■ Mt X (C*)^ ~ M^'^ xT^ with A^f = Higgs scalars in the fundamental representation. In the case 
. of the Abelian-Higgs model {N-p = Nq = 1), the intersecting vortex sheets can be beautifully 
understood in a mathematical framework of amoeba and tropical geometry, and we propose 
a dictionary relating solitons and gauge theory to amoeba and tropical geometry. A projective 
shape of vortex sheets is described by the amoeba. Vortex charge density is uniformly distributed 
among vortex sheets, and negative contribution to instanton charge density is understood as the 
complex Monge- Ampere measure with respect to a plurisubharmonic function on (C*)^. The 
Wilson loops in are related with derivatives of the Ronkin function. The general form of the 
Kahler potential and the asymptotic metric of the moduli space of a vortex loop are obtained as 
a by-product. Our discussion works generally in non-Abelian gauge theories, which suggests a 
non-Abelian generalization of the amoeba and tropical geometry. 
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1 Introduction 



The study of topological solitons is intimately connected with the development of mathemat- 
ics. Originally solitons are found as solutions to non-linear (field) equations in various physical 
systems, but many mathematical concepts have been developed at the same time in order to un- 
derstand properties and integrability of the soliton equations. Indeed, many mathematical tools, 
such as representation theory of differential operator algebra and infinite dimensional Grassman- 
nian, are required to investigate the solitons. The language of two-dimensional conformal field 
theory is also useful in studies of a class of the soliton systems, such as KdV, Toda, KP equations. 
The investigation of solitons in gauge theory, such as instantons and monopoles, gives not only 
non-perturbative informations about the field theory, but also stimulates new developments of 
differential geometry. 

The purpose of the present article is to study topological solitons using novel mathematical 
objects known as amoeba and tropical geometry. We find a one-to-one correspondence between 
amoeba/tropical geometry and solitons in Yang-Mills-Higgs theory. Yang-Mills gauge fields cou- 
pled to the Higgs fields naturally appear as a bosonic part in supersymmetric Yang-Mills theory 
with eight supercharges. When sufficient number of the Higgs fields get vacuum expectation 
values (vevs), the theory is in the Higgs phase with completely broken gauge symmetry. Typical 
solitons in the Higgs phase known for a long time are vortices in the Abelian U{1) gauge the- 
ory coupled to a single complex Higgs field (the Abehan-Higgs model) These vortices have 
been recently extended to the non-Abelian case, vortices in completely broken non- Abelian U{N) 
gauge symmetry [2]-[T2]. The other fundamental solitons found relatively recently in Yang-Mills- 
Higgs theory (or corresponding nonlinear sigma models) are domain walls or kinks [6], [T3] |16]. 
Since both vortices and domain walls preserve half of the supercharges when embedded into su- 
persymmetric theories, they are called 1/2 BPS solitons. Composite solitons in the Higgs phase 
have recently been studied extensively, especially in supersymmetric gauge theories with eight 
supercharges [171 UHl [19]. Since the magnetic field has to vanish in the Higgs phase, magnetic 
monopoles are confined by vortices (confined monopoles) [2n]^[23]. Although isolated instan- 
tons shrink to points in the Higgs phase, they can lie inside a vortex core (trapped instantons) 
[22l [23] . Vortex-strings can end on a domain wall [24j or stretch between domain walls [25] . 
When a (composite) soliton configuration breaks n directions of translational symmetry, it is 
defined to have n codimensions. The composite soliton with the lowest co-dimension is a web 
or a network of domain walls, whose moduli space and dynamics have been worked out recently 
[26]- [30]. All of these composite solitons preserve 1/4 of supercharges if we realize the Yang- 
Mills-Higgs theory as a supersymmetric Yang-Mills theory, and are called 1/4 BPS states. It has 
been found that all these 1/4 BPS composite sohtons are related by the (Scherk-Schwarz twisted) 
dimensional reduction, starting from the instanton- vortex system [23]. Therefore understanding 
the instanton-vortex system is of primary importance to study BPS solitons in supersymmetric 
Yang-Mills-Higgs theories with eight supercharges]^ Unfortunately, generic configurations of in- 

^ It is worth pointing out that there exist the other series of 1/4 BPS systems. This contains a triple intersection 
of vortices where a set of two vortices has one common codimension [31j , contrary to the instanton-vortex system 
|23| where vortices have no (or two) common codimensions. The former preserves (1,1) SUSY and the latter (2,0) 
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stantons and vortex sheets as co- dimension four solitons have not been worked out, apart from 
a trapped instanton as a lump on a (uncurved) vortex plane, or an intersection point of two 
orthogonal (uncurved) vortex planes [23|. It is important to characterize generic configurations 
of instantons and vortex sheets as co-dimension four solitons in a precise and transparent man- 
ner. We can call these solitons as webs of vortices. We show that these soliton webs are nicely 
described in terms of the amoeba and the tropical geometry. 

Interestingly, amoeba and tropical geometry have already appeared in physics literature, in 
the context of topological strings. Topological string amplitude (or its building block, topological 
vertex) can be described by means of a melting crystal picture [3^ [35] , which was originally 
introduced as a statistical model in mathematical physics. Shape of the melting crystal in 
the thermodynamic limit is well interpreted in terms of an amoeba, which is a logarithmic 
projection of a smooth Riemann surface. On the other hand, an emaciated body of the amoeba 
corresponding to zero temperature limit can be captured by tropical geometry, where we can 
obtain some properties of the Riemann surface from a skeleton of the amoeba. In the context of 
superstring theory, the amoeba and tropical geometry appear as, respectively, quantum ((/^ ^ 1) 
and classical [qs — >■ 0) shape of intersecting five-branes (five-brane junctions), which are dual to 
a suitable Calabi-Yau geometry where the topological string is defined. Other examples appear 
in an instanton and BPS state counting problem of the supersymmetric (quiver) gauge theories 
[36| [371 [38] . since these systems are closely related to the topological string amplitudes and 
realized in terms of the five-brane web. 

The five-brane system is useful to understand the relationship between the gauge and string 
theory, and the amoeba and tropical geometry, which is the main subject of this article. Our 
study is inspired by these successful applications of the amoeba and tropical geometry to physics. 
We find a similarity between the five-brane web and the 1/4 BPS composite solitons of vortex 
sheets and instantons. Some of properties of the webs of vortices are still unclear. In particular, 
the instanton charge appears at the intersection point of the vortices, but we have not obtained 
a tool to see a distribution of the extra instanton charge on the vortex web. The five-brane web 
in the superstring theory and vortex web possess common properties. Indeed, toric diagrams 
and geometry play an important role in both sides, and their moduli spaces are described by 
similar quotient spaces and the moduli space of the vortex web should be included in that of 
the five-brane web. So we expect that the amoeba and tropical geometry, which is important to 
understand dynamics of the five-brane web, are also useful to analyze the web of the solitons in 
the Yang-Mills-Higgs system. 

In this paper, we study the most generic configurations of 1/4 BPS solitons of instantons 
and vortices in the Higgs phase of the five-dimensional A/" = 1 supersymmetric U{Nc) gauge 
theory (with eight supercharges) on x (C*)^ ~ M."^'^ x with Np = Nq Higgs scalars in the 
fundamental representation, by using the moduli matrix formalism [18]. Torus T"^ = x 
allows us to obtain other 1/4 BPS solitons readily through dimensional reduction. We show 
that vortex sheets are defined by zeros of a Laurent polynomial of two complex coordinates of 

SUSY in terms of two-dimensional supersymmetry [32]. Different sets of 1/4 BPS equations in these two series 
are obtained as dimensional reductions of the unique set of the 1/8 BPS equations [32l [33] , 
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(C*)^, and instanton positions are given by common zeros with another polynomial. We find that 
the above expectation of the importance of amoeba and tropical geometry to be correct. We 
describe physical quantities of the intersecting solitons (soliton web) in terms of the mathematical 
language of amoeba and tropical geometry. We also see the important objects in the amoeba and 
tropical geometry, such as the logarithmic mapping, the Ronkin function and the Monge- Ampere 
measure, also have essential meanings on the soliton side. We describe properties of the soliton 
web in terms of these mathematical objects. We find that the moduli matrix approach is very 
useful in the translation between physical and mathematical languages. 

The organization of the present paper is as follows. In Sec. [21 we review the 1/4 BPS equa- 
tions for vortices and instantons in the Higgs phase of five- dimensional A/" = 1 supersymmetric 
U{Nc) gauge theory on Rt x (C*)^ ~ M^'^ x with Np = Nq Higgs scalars in the fundamental 
representation, using the moduli matrix formalism. We mainly consider vortices in the overall 
U{1) gauge theory except in Sec.[5l In Sec.[3l we consider simpler case of vortex on cylinder 
R X S'"'^ ~ C*. In Sec.Hl we study the most general situation of vortex sheets on (C*)^. Sec. 14. II 
gives amoeba corresponding to vortex sheets. Sec. 142] relates it to the tropical geometry and gives 
an example. Sec. 14. 31 describes a general formula to compute the topological charges. Sec. 14. 41 
gives the metric of the moduli space of a vortex loop. We also obtained a new general formula 
for the Kahler metric which is valid for arbitrary values of Nq and Np. In SecO we exam- 
ine the instanton number in non-Abelian gauge theory. Sec. 15. II reviews the instanton number 
for instantons trapped inside a non-Abelian vortex plane [23]. Sec. 15. 21 discusses more general 
configurations of the instantons trapped inside a non-Abehan vortex web. Sec. [6] is devoted to 
conclusion and discussion. 



2 Vortices and Instantons 

In this section, we first review the construction of vortex solutions in U{Nc) gauge theory in 
(4-|-l)-dimensional spacetime x (C* with Np Higgs fields in the fundamental 

representation. By introducing additional A^f Higgs fields in the fundamental representation, 
this theory can also be regarded as the bosonic part of a five- dimensional A/" = 1 supersym- 
metric U{Nc) gauge theory with Np hypermultiplets in the fundamental representation, but the 
fermionic part (and another set of Np Higgs scalars) is irrelevant and is omitted in the following 
discussion. We introduce (xi, yi, X2, 1/2) and zi = xi + iyi, Z2 = X2 + iy2 as real and complex 
coordinates of (C*)^, respectively. 

The Lagrangian of the theory takes the form 



£ = Tr 



2g'- 



(2.1) 



where the Higgs fields are expressed as an A^c X Np matrix H"-^ {r = 1, ■ ■ ■ , Nc, A = l,--- ,Np). 
The covariant derivative is defined by P^-fT = d^H + iW^H and the field strength by F^^ = 
—i[D^, Vy] = dfj_W^ — d^W^ + i[W^, Wy\. The constants g and c are the gauge coupling constant 
and the Fayet-Iliopoulos (FI) parameter, respectively. At the vacua (the minima of the potential) 
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of this theory, the Higgs fields H get vev and U{N) gauge symmetry is completely broken. Namely 
the theory has only the Higgs branch due to the nonzero FI term. The moduli space of the vacua 
is given by a complex Grassmannian 

= SU{Nc)xsf{Nj-Nc)xU{iy ^^'^^ 

Considering a static gauge configuration, we find that there is a lower bound of the energy 
[221123] 

E> - ^ [ Ti(F AF) -c [ TiF ALU = ^I + 2ttcV, (2.3) 
J J 9^ 

where the two form uj = ^{dzi A dzi + dz2 A dz2) is the Kahler form on (C*)^, and we have defined 
the total instanton charge / as an integral of the instanton charge density X, and the vortex 
charge V a,s a. divergent integral of the vortex charge density V 

I ^ Jl ^ -±-jTriFAF) = Jch, (2.4) 

V = [v = f TiFAuj = ! CiAu. (2.5) 



27r , 

The lower bound Eq. (12. 3 P is saturated if the following BPS equations ^ 

F,,,, = 0, V,^H = 0, -2i{F,,,, + = ^(m^ - cl^J, (2.6) 

are satisfied. When FI parameter c is sent to zero, the Higgs field H vanishes and these equations 
reduce to the anti-self-dual equations for Yang-Mills instantons, whereas when we neglect the 
2:2- (or zi-) dependence of the fields they reduce to simple vortex equations for vortices on the 
2;i-(or Z2-)plane. These vortices are two- co dimensional surfaces in the four dimensional space 
(C*)^. Therefore the BPS equations (12.61) contain at least instantons and intersecting vortex 
sheets. As we will see below, these equations describe webs of vortex sheets in general. The 
equations (12. 6p were earlier found for those on arbitrary Kahler manifold [39j and were simply 
called "vortex equations" although they contain instantons also|^ It has been shown in [23] that 
solutions to the BPS equations ^Mj on (C*)^ (or C^) preserve a quarter of supercharges in the 
supersymmetric gauge theory with eight supercharges. So the configuration of the solution is 
called a 1/4 BPS state in this sense. The energy of the BPS configuration is determined by the 
topological charges (12.41) and (12. 5p . 

The vortex charge V can be evaluated from the sum of the area of each vortex sheet, as we 
will see in Sec. HI The total instanton charge I can be decomposed into the intersection charge 



^ In Ref. |39| the vortex equations are defined on arbitrary Kahler manifold M of complex dimension n (with 
n = 2 not necessary). There the bound is given by J^^ Tr {F A F) A and J^^ Tri^ A oj""^ with the Kahler 



2-form uj, instead of the charges (|2.4[) and (|2.5p . Furthermore at least in the case of A^'c ~ Np ^ 1 these generalized 
equations can be obtained as equivariant dimensional reduction of the Donaldson-Uhlenbeck-Yau equations on 
M X with a monopole configuration on S*^ [40j . 
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/intersection and the instaiitoii number Jinstanton as 



^ -^intersection ~l~ -^instanton; (^■'^) 

-^intersection = ^ ^^intersection = JttFATtF = - J CiACi, (2-8) 

-^instanton = / -^instanton = / ^2- (2.9) 



The intersection charge has negative contribution to the energy of the BPS configuration, which 
can be regarded as binding energy of intersecting vortex sheets. On the other hand the instanton 
number (not to be confused with total instanton charge) counts to the number of usual (particle- 
like) instantons and has positive contribution to the energy. 

Let us solve the BPS equations (12.61) . The first BPS equation -Fk^sa = "^[^21)^^22] = in 
(12.61) is equivalent to an integrability conditioijfl for the differential operators Vg^ , which states 
that there exists an A^'c x Nq matrix-valued functioii^ S{zi, zi) G U{Nq)'^ = GL{Nq, C) such that 



W,^ = -iS"'d,.^S. (2.10) 



Defining an A^^c x N-p matrix 



= SH, (2.11) 



the second equation in (12. 6p reduces to 



d,,Ho = 0. (2.12) 

This means that the elements of the matrix Hq should be holomorphicEl with respect to the 
complex coordinates Zi. The matrix- valued quantity S is determined from the last equation in 
Eq. (12. 6p . which can be rewritten in terms of an A^^c x A^c positive definite Hermitian matrix!^ 

Q = SS^ (2.13) 

into 

2 

d,,{nd,,n~') + d,,{nd,,n-') = {i^^ - no^-') , (2.14) 

where we have defined 

no = \hoHI (2.15) 



^ If we identity the Higgs field H"^^ as a set of A'f sections of a rank Nq vector bundle E on the base Kahler 
manifold, then the first BPS equation i^ziZ2 = is equivalent to the condition of the existence of the holomorphic 
frame {e^} (i = 1, • • • , iVc), which satisfy V^ei = 0. 

^ The complexified gauge transformation iSV can be interpreted as the change of basis from the unitary frame 
(orthonormal frame) {e^} (r = 1, • • • , Nq) to the holomorphic frame {ci}. 

^ In other words, Hq"^ is a set of A^f holomorphic sections of the holomorphic vector bundle E. 

^ This matrix f2 can be interpreted as the inverse of the Hermitian metric in terms of the holomorphic frame 
{ci} which is the identity matrix in terms of the unitary frame {e^}. 
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We call the equation f l2.14p the "master equation" of the instanton-vortex systemjll 

Using these redefined fields, we can solve the BPS equations by the following procedure. 
Take an arbitrary holomorphic matrix Hq{z) and solve Eq. fl2.14p in terms of Q, then we can 
determine S up to U{Nc) gauge transformation S SW and physical fields can be obtained 
via the relations 

W,^ = -iS-'d,^S, H = S-'Ho. (2.16) 

The equations (12.121) and (12.141) have a "gauge symmetry", which we call "V^-transformation" , 
defined by 

{Ho,S) ^ {VHo,VS), V{z)eGL{NcX)- (2.17) 

Note that the physical fields and H are invariant under the V^-transformation, so this defines 
an equivalence relation called the "V^-equivalence" Jf| Assuming that there exists a unique solution 
of Eq. (12.141) for a given -ffo(2;)|^ we find that there exists a one-to-one correspondence between 
the equivalence class Hq ~ VHq and a point on the moduli space of the BPS configurations. In 
this sense, we call Hq{z) a "moduli matrix" and the parameters contained in Hq are identified 
with the moduli parameters of the BPS configurations. 

Now let us consider the case Nq = Np = N which is often called a local theory. If the 
determinant of the Higgs fields H vanishes in some regions, the broken gauge symmetry is partially 
restored in those regions. As dictated by the Meissner effect in the Higgs phase, this gauge 
symmetry restoration occurs where the magnetic fiux penetrates and a vorticity arises around the 
zero of the Higgs field. Therefore the vanishing determinant det H = defines a two-dimensional 
surface of vortex positions in the four-dimensional space: it is equivalently given b}0 

detHo{zi,Z2) = 0. (2.18) 



In closing this section, we comment on the vortex solutions of the (2-|-l)-dimensional U{Nc) 
gauge theory on x C* with Np massless Higgs fields in the fundamental representation, since 
we study this case in the next section. Historically, the non-Abelian vortices were first found [21 [3| 
on C in the color-fiavor locked phase of U{Nc) gauge theory with Np massless Higgs fields in the 
fundamental representation. After their discovery the non-Abelian vortices have been extensively 
studied by many authors In particular, the non-Abelian vortices on a cylinder C* have been 
studied in [TJ ^J]. In the moduli matrix formalism, the discussion of the C* case is completely 

When the Fl-parameter c goes to zero, the Higgs phase no longer exists. In this case, the RHS of Eq. (12.141) 
vanishes, and Eq. (|2.14p becomes the so-called Yang's equation 41] for usual instantons not accompanied by 
vortices. 

* This equivalence relation originates from the redundancy of the holomorphic frame {ii}. 

^ This assumption is correct at least when the base space is compact and Kahler, since the uniqueness and 
existence of solutions to the BPS equations (|2.6p were rigorously proved in terms of the Hitchin-Kobayashi 
correspondence 

^'^ The determinant detiJo can be regarded as the holomorphic section of the determinant line bundle A^'^E, 
and the vortex sheet corresponds to the effective divisor associated with the holomorphic section det Hq . 



6 



parallel to the (C*)^ case and the necessary formulae are simply obtained by neglecting the Z2, Z2 
dependence. For example, the BPS bound (12. 3p is reduced to 

E>-cl^.TrF.„ (2,19) 

and the master equation fl2.14p for Q{zi) is given by ^ 

2 

d,,{nd,,Q-') = (1^, - Qo^-') . (2.20) 

For the Abelian-Higgs model, Nc = Np = 1, this equation reduces to the so-called Taubes's 
equation [42] after some redefinition. 



3 Vortices on a Cylinder C* 

In this section, we first consider the simpler case of vortices on a cylinder R x ^ C* be- 
fore discussing the intersecting vortices in four dimensions. The vortices are BPS solutions of 
the (2+l)-dimensional U{Nc) gauge theory with Np massless Higgs fields in the fundamental 
representation, with one spatial direction compactified [71 Hi]. In the supersymmetric system, 
the vortices preserve a half of the supercharges. Although our real interest is in complex two- 
dimensional case (M^ x ~ (C*)^), one-dimensional case is simpler and useful to understand 
the discussion in the next section. 

Let (x, y) and z = x + iy he real and complex coordinates of R x 5*^ ~ C*, respectively. The 
coordinate of 5*^ has a period 27ri?, namely y ^ y + 2ttR. We here concentrate on the case of 
Nc = Np = N. The moduli matrix formalism works as well in this (2+l)-dimensional case, and 
the BPS solutions are parametrized by the moduli matrix Hq{z). Since the moduli matrix should 
satisfy the periodic boundary condition Hq{z + 2mR) = Ho{z), the determinant of the moduli 
matrix can be expanded as a Fourier series 

detHo = ^ant""'!^. (3.1) 

Introducing a new coordinate u = e^^^, this can be rewritten as 

detHo = P{u) = ^a„M". (3.2) 

The positions of the vortices are determined by zeros of this Laurent polynomial. By performing 
an appropriate ^^-transformation Hq{z) V{z)Ho{z), detHo reduces to 

k k 

detHo = Hie'/'' - e'^/'') = l[{u-u,), = e^»/«, (3.3) 

i=l i=l 

with k denoting the number of the vortices. The solution Q of the master equation (12.201) for the 
vortices approaches to Qq = -HoHq in the strong gauge coupling limit g ^ oo. In this limit the 
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configuration of tlie magnetic flux of the overall U{1) becomes singular such as 

k 

TrF = i(9(91ogdetfi idd log\ det HqI"^ = -27t^ 6^ {z - Zi)dx A dy, (3.4) 

i=l 

with d = dzdzi d = dzdz- This reflects the fact that the size of the vortices is proportional 
to / = \lg\fc and becomes zero in the strong coupling limi10. Using the configuration in the 
infinite coupling limit, we find that the topological charge is given by the number of zeros of the 
polynomial Piu) 



1 

2^ 



d^xTTF^y = k. (3.5) 



If we dimensionally reduce the theory on S^, then the vortex can be viewed as a domain wall 
in (l+l)-dimensional theory [14]. These field theoretical BPS solitons are realized by kinky D- 
brane configurations in superstring theory [T51 [TB] , and the relation between vortices and domain 
walls is understood via T-duality [7j. In field theory language, the profile of the kink solution of 
the domain wall is described by a logarithm of a Wilson line along 



Pexp ^ J dyWy 



(3.6) 



S(a;) can also be viewed as the adjoint scalar in the T-dual (dimensionally reduced) theory. See 
Fig. Wljo) for an example of a Tr S(a;) plot. 

It is convenient to define a function Np{x) associated with the Laurent polynomial P{u) = 
det Ho{z) by 



f^""^ dy 1 
Np{x) = lim / — — -logdetn 
g^oo Jq ZnK I 

"2-^ dy 



2-nR 



log I det i^ol 



1 ^ 

-^^^{xd{x - x.^ \ Xid{xi- x)^. (3.7) 



In the final line, we have used the Jensen's formula |43j^n . As we will see in the next section, 
this piece-wise linear function is the "Ronkin function" in one dimension. By using this function 



^^The appearance of length scale I = \/g^ is understood from the master equation (|2.20p . In (|2.20p . the 
parameters g and c appear only in the combined form g^c, and thus / is the only length scale made from g^c. 

^^The classical Jensen formula states that for arbitrary holomorphic function f{x) with zeros at at {i — 
1,2,... k), we have 

^ /'^og|/(re^'^)|rf0 = log 1/(0)1 + ^ log l-l, (3.8) 
277 7o ^ a, 

where we have chosen indices i such that ai < r for i = 1,2,... Nr and ai > r otherwise. 
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Np, the trace of the adjoint scalar in the infinite gauge couphng hmit can be written by a step 
function 



hm TrS = d^Np{x) 

g~*oo 



i=l 



Xi = -Rlog|Mj| = ReZi. 



(3.9) 



Note that in the strong gauge couphng hmit, the smooth kinky profile reduces to step- wise 
profilJ^ as shown in Fig. [1] (c). 




2ttR 





11 R i; 


IjR 






1/7? 








xi 


Xi 



(a) energy density 



(b) Tr t 



(c) lim Tr S = d^Np{x) 

g^oo 



Fig. 1: (a) represents the energy density of two vortices. The energy is localized around the 
center of the vortices, (b) shows the profile of a kink solution, or equivalently Tr S as defined in 
(13.61) . In the strong gauge coupling limit, the profile reduces to step- wise shape as shown in (c). 



We can also express topological charges in terms of Np and its derivative S. The BPS bound 
(12.191) is now rewritten as 



E>-c J d^xTiF^y = cj dxd^TiJ:, c = 2ttRc, (3.10) 

where we have used / dy dyTi Wx = and Tr S = f dy Tr Wy. In the context of the domain 
wall, the quantity c / (ixS^TrS provides the sum of the charges (masses) of the domain walls 
|14j . Interestingly, the energy of the BPS configuration can be determined only from Tr E, namely 
the zero mode of the overall U{1) gauge field TiWy. This is because the topological charge is 
determined only from the boundary condition and all the massive KK-modes vanish at spatial 
infinities x ±oo in the BPS configurations. 

We finally discuss the Kahler metric on the moduli space of BPS vortices. The moduli space 
of the vortices is a Kahler manifold parametrized by the moduli parameters (pi, (pi- The Kahler 
metric of the moduli space is directly calculated from the solution of the master equation as 
follows [S] (see also [231 [IS]) 



ICfj = Tr 



J d'^x }Cij{z, z, R,l), (3-11) 
didjlogn + Al^(dg{ndin-^)dj{nd,n-^) - ds{nd^n~^)dj{ndin'^))] , (3.12) 



This limit is different from the one taken in [7\ where the profile is not step function but has a constant slope 
in the interval of a vortex size. 
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where di = d/d(f)i, dj = d/d(j)j are derivatives with respect to the moduh parameters and 
/ = \j is the length scale of the vortex core. Note that the matrix- valued function f2 depends 
on the parameters g and c only through /. 

For concreteness, let us consider fc-vortex configurations in the Abelian-Higgs model {Nq = 
Np = 1). In this case the 2/c-dimensional moduli space is parameterized by the positions of 
vortices Zi {i = 1, ■ ■ ■ , k). For well-separated vortices \zi — Zj\ R, I, the asymptotic metric is 
obtained by taking the limit I —>■ 0, R —>■ 0. In the small vortex limit / — > the Kahler metric 
becomes 

K,^ ^ c / d''x—--\og\Ho\\ (3.13) 



dzidzj 

Therefore the Kahler potential which determines the metric by a relation Kf^ = q^q- , can 
be written by 

K ^Attc j dx (^Fp{x, Zi, Zi) - /(x, Zi) - f{x, Zi)^ , (3.14) 

where Fp{x) = lim/j^o RNp{x) and f{zi) is a holomorphic function which is required to make the 
Kahler potential finite. This Fp is a one-dimensional "tropical polynomial" that we will extend 
to the two-dimensional case in the next section. Since the asymptotic forms of the function Fp[x) 
are given by 

17 f \ j Xi + --- + Xk, X ^ -oo 

^-^"^ = 1 kx, x-oo ' ^'-''^ 

a possible choice of the function f{x, Zi) is 

/(x, z,) = ''^"^^'' 9i-x) + ^9{x). (3.16) 

Then the asymptotic Kahler potential can be evaluated as the area of the shaded region in Fig. [21 
which is given by 

27^0^2^', = ^f2('^ + -z.r. (3.17) 
1=1 1=1 

Using this Kahler potential, the Kahler metric is given by K^^j = didjK = T[c5fj. So the effective 
Lagrangian which describes the dynamics of well-separated vortices becomes 

k 

LeS = KfjZiZj = TTC^ \Zif. (3.18) 

i=l 

This shows that the well-separated vortices behave as undistinguished free particles with the 
mass 27rc. 

In the next section, we move on to the case of vortex on (C*)^. Although the story is more 
complicated and general, we will encounter similar structures to those described in this section. 
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Fig. 2: The asymptotic Kiihler potential can be evaluated as the area of the region surrounded 
by Fp{x) and f{x,Zi) + f{x,Zi) (shaded regions). The area of the meshed region gives the 
contribution from the center of mass modulus ^{z^ + z^)'^, Zc = {zi + Z2 + ■ ■ ■ + Zk)/k. 



4 Webs of Vortex Sheets on (C*^^ 
4.1 Vortex Sheets and Amoeba 

Let us consider the vortex- instanton system on (C*)^ ~ x T^. As before, we will use 
(xi, yi, X2, 112) and zi = xi + iyi, Z2 = X2 + iy2 as real and complex coordinates of (C*)^, respec- 
tively. The coordinates of are identified with periods (27r_Ri, 271 R2), namely yi ^ yi + 27r_Rj. 
In this case, the determinant of the moduli matrix det Hq, which defines the vortex sheets, is 
written in the form of the Fourier series 

deti/o(^l,^2) = «nl,n2e^'^ + ^'^ (4.1) 

(ni,n2)eZ2 

If we define new cylindrical coordinates (mi,M2) on (C*)^ by Ui = e-^*, det i^o is now written by 
a Laurent polynomial 

P{uuU2) = det Ho = Yl «m,n2wrwr- (4.2) 

(ni,n2)GZ2 

The positions of the vortices are described by zeros of P{ui,U2) similarly to those on C* in the 
previous section, but the vortices form a two-dimensional sheet (surface) in (C*)^ in the present 
case. 

We define the "Newton polytope" A(P) C of a Laurent polynomial P{ui,U2) by 

A(P) = conv. hull |(ni,n2) G Un^n^ ^ 0} . (4.3) 

Conversely, P{ui,U2) is called the Newton polynomial of A, when its Newton polytope A is 
convex. In the discussion of domain wall webs, the Newton polytope A(P) was called the 
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"grid diagram" [26]-[30]. When we say "Newton polynomial", the coefficients ani,n2 in (14. 2p are 
arbitrary parameters. Namely, a„j^„2 regarded as moduli parameters of the vortices. 

Analogous to the case of the vortices on the cyhnder discussed in the the previous section, a 
web of vortices on (C*)^ is now dimensionally reduced to a web of domain walls on [26] |30]. 
In order to see the connection better, we define "amoeba" of P b}0 

= |(i?llog|Mi|, i?2log|M2|) G I P(ui,M2) =0}. (4.4) 

Note here that Ri log = Xi and R2 log \u2\ = X2- This is a projection of the shape of vortex 
sheet onto two non-compact directions. See Fig. [3] for an example of amoeba. From this example. 




(a) Newton polytope (b) amoeba 



Fig. 3: An example of amoeba; P{ui,U2) = ao,o + ai.o^i + cl2,oUi + ct3,o^i + cto,i^2 + ai,iMi^i2 + 

a2,lM?M2 + a.3,lM?M2 + 00,2^1 + CH,2MlM2 + «2,2^i?M2- 

we can see that the amoeba has several asymptotic regions extending to infinity, which are called 
"tenticles" in literature. In physics terms, each of these tenticles represents a semi-infinite cylinder 
of the vortex, and, for generic coefficients, corresponds to the normals to the Newton polytope 
(see Fig. [3])o We also have (again, for generic values of moduli) holes for each internal lattice 
point of the Newton polytope. We learn from these facts that amoeba is a projection of generic 
webs of the vortices. 

The notion of amoeba dates back to It was originally studied in the context of mon- 
odromy of the so-called GKZ-hypergeometric (originally called ^-hypegeometric) functions. It is 
also intimately connected with real algebraic geometry (Hilbert's 16th problem). Furthermore, it 
plays an important role in the discussion of the tropical geometry, as we will see. What we have 
found is that amoeba also appears quite naturally in the discussion of the webs of the vortices. 
See also the work [38], which discusses amoeba in the context of instanton counting. 

To see the relation between the webs of the vortex sheets and the webs of the domain walls. 

Amoeba can be defined for (C*)" with arbitrary integer n, but we only use the case of n = 2. 
^^When moduh parameters are chosen to be special value, several tenticles of amoeba can merge into one. In 
this case, the multiplicity of the spires is considered to be greater than one. 
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it is convenient (just as in the previous section) to define Si(xi,X2) by 



dy2 

2ttRi J 2mR2 



log 



P exp 



dyiW, 



yi 



(4.5) 



and similarly for Tj2{xi,X2) by interchanging the subscript 1 with 2. These adjoint scalar fields 
i]j(xi, X2) on are interpreted as the zero modes of the gauge fields in the Kaluza-Klein decom- 
position and exhibit the two-dimensional kink profiles. The trace of these adjoint scalar fields 
can be simply written as 



Tr 



Si(a;i,X2) 



1 



1 



2'kRi 27ri?2 
1 d 
8it'^RiR2 dxi Jj-2 



(fylog detn, 



(4.6) 



where we have used Eq. (12.101) and fl = SS'^ . The matrix- valued function fl is simplified in the 
strong gauge coupling limit to (recall master equation (12.141) ) 



lim 



-HoHl 



(4.7) 



Although the vortex sheet becomes thin and singular in the strong coupling limit, the matrix Qq 
still has important physical informations on vortices and instantons. The traces of the adjoint 
scalar fields Eq. (14.61) in the strong coupling limit are given by 



lim Tr 

g^oo 



ti{xi,X2) 



^ AT / X 

— A/p(Xi,X2 
OXj 



(4i 



where Np{xi,X2) is nothing but the "Ronkin function" [45j in two dimensions 

1 1 



Np(xi,X2) 



27iRi 2nR2 Jt'2 
1 



d y\og\ detHo{zi, Z2) 



dui du2 , I „ , . I 

A log|P(Mi,U2)|, 



(4.9) 



(27ri)2 

defined from the Laurent polynomial P{ui,U2) = det Ho{zi, Z2). 

The Ronkin function has several interesting properties. First of all, it is convex [15]. Second, 
the derivatives of the Ronkin function, Tr Si and Tr S2, take constant values in each complement 
of the amoeba0 and those constant values (multiphed by R1R2) are given by the lattice points 
in the Newton polytope of P [16]. More generally, i?ii?2(Tr Ei, Tr E2) as a function defined on 

(including points on the amoeba) take values within the Newton polytope A(P) of P. 



4.2 Relation with Tropical Geometry 

Now one difference arises from the previous section. In the complex one- dimensional case dis- 
cussed in the previous section, the Ronkin function Np is piece-wise linear when we take the 

^^In [35] and many other literature, these constant values are called the orders of the complement of the amoeba. 
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thin wall limit / = 1/ g^/c — > 0. In two-dimensional case, however, the Ronkin function and its 
derivative are smooth even when the gauge coupling goes to infinity. 

We can still consider another limit in which the derivative of the Ronkin function becomes 
discontinuous. The limit is i?i = i?2 = -R with fixed 

rm,n2 = ^log|a„i,„2|. (4.10) 

This limit corresponds to dimensionally reducing the theory to (2+l)-dimensions, neglecting all 
KK modes. In this limit, the amoeba degenerates into a set of lines ("spines"), which is called 
"tropical variety" in the tropical geometry literature. Physically speaking, the vortices reduce 
to the domain walls by dimensional reduction, and the tropical variety signifies the location of 
the domain walls. At the same time, the Ronkin function becomes a piece-wise linear function 
Fp{xi,X2) defined by 

Fp{xi,X2) = lim i? log |P(mi,-U2)| = max (t^i^i + 7^2X2 + r„^,„2) , (4.11) 

where V{Q) is a set of the vertices associated with the Newton polytope Q, and r„j in Eg. (14. 101) 
are constants determined from constants a„j^„2 in (14. 2p . This function Fp coincides with the 
Ronkin function Np on the complement of the amoeba (recall Np is linear on each complement). 

If we compare (14.111) with (14. 2p . we notice that the sum and products in the polynomial 
Yl n2'Wi^^2^ are replaced by a maximum function max (niXi + ^2X2 + of the lin- 

(ni,n2)ey(Q) 

earized functions. The formal reasoning is given as follows. If we define xi = exp{xi/ R), 
X2 = exp{x2/ R), Xi + X2 = exp{x2,/ R), and X1X2 = exp^x^/ R), then we find 

Xs = max(a;i, 0:2), X4 = Xi + X2 (4-12) 

in the R ^ 00 limit. Hence in the tropical limit, the ring (R, -|-, x) is replaced by an idempotent 
semirinj^ (M, ©, ®), with a tropical addition © and a tropical multiplication © given by 

xi © a;2 = max(xi, X2), xi ® X2 = xi + X2, (4-13) 

respectively. The semiring (M, ©, ©) is sometimes called the tropical semiring or the max-plus 
algebra. 

The operation replacing the addition and multiplication with the tropical addition and tropical 
multiplication is also called dequantization or ultradiscretization. It appears in a discretization of 
integrable soliton equations such as KdV, Toda and KP hierarchies and also in cellular automata. 
These integrable soliton systems seem to be completely different from the vortex-instanton system 
we are considering, but it is interesting that the same structure plays important roles in many 
integrable systems. 

We have mentioned about the tropical limit and tropical semiring, but then what is the 
corresponding geometry? In usual algebraic geometry, we consider geometry corresponding to 

"'^^A semiring is an algebraic structure similar to a ring, but without the requirement that each element must 
have an additive inverse. 
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commutative ring. In contrast, the geometry corresponding to tropical semiring (M, ©, (g)) is 
called tro,,cal (^,ebra,c) geornetryQ we can formulate and prove "tropical analogue" of many 
theorems in usual algebraic geometry, such as the Riemann-Roch theorem and the Bezout's 
theorem. Although the study of idempotent semirings in applied mathematics (such as control 
theory and optimization) has a long history [50], the study of corresponding geometry is relatively 
new and it is still an active area of research (see [511 H?! El])- The tropical geometry has 
now diverse applications, ranging from string networks [53], enumeration of curves [51], mirror 
symmetry [55] and even computational biology [56] . 

We can consider a tropical version of algebraic variety, namely tropical variety. In the liter- 
ature, it is often defined as a non-Archimedian amoeba, but for our applications, it suffices to 
define it as the set of points where the piece-wise linear function Fp{xi,X2) ("tropical polyno- 
mial") is not differentiable. This is nothing but the skeleton (spine) of the amoeba in the limit 
i? — i> 0, and its physical meaning is the position of the domain walls, namely the position of the 
step-wise kinks appearing in the profiles of TrSj(xi,X2) {i = 1,2). An example of the tropical 
varieties are shown in Fig. HI As shown there, tropical varieties (in the situation we want to 
consider) are obtained from triangulation of the Newton polytope ([51], Proposition 3.5). In 
this sense, this is similar to the so-called (p, g)-web or web diagram in [57]. We will make more 
comments on this analogy in the last section devoted to the discussion. 

(a) amoeba 

Fig. 4: An example of the amoeba and corresponding tropical variety. 

Example 

Let us here give a simple and concrete example for later discussions. An example of the Newton 
polytope is given in Fig. [5] (a) and we set Ri = R2 = 1 for simplicity in the following. Then the 
corresponding Newton polynomial is given by 

P(ui, U2)=ui + U2 + 1 = + e'' + 1, (4.14) 

and its amoeba and corresponding variety is shown in Fig. [S] 

"'^^ According to [47j, the name "tropical" was coined by a French mathematician Jean-Eric Pin [48 , in honor of 
their Brazihan colleague Imere Simon [15] . 
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(b) tropical variety 
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(a) Newton polytope 



(b) amoeba 



(c) tropical variety 



Fig. 5: The Newton polytope (a) amoeba (b) and tropical variety (c) for the Laurent polynomial 

P(ui, U2) = Ui + M2 + 1. 



The derivatives Tr Si(xi,X2) of the Ronkin function associated with P are computed to be 



lim Tr 



Si(xi,X2) 



1 COS 

TT 



-1 







for xi < log |e^'^ — 1| 

for log \e^^ — 1| < Xi < log {e^"^ + 1| 

for xi > log \ e^'^ + 1| 



lim Tr 



S2(Xi,X2) 



1 





1 , f e2^2 _ gSxi _ ^ 
1 cos ' 



TT 



for X2 < log I e^'^ — 1 1 

for log le""^ - 1| < a;2 < log le""^ + 1| 

for X2 > log I e""'^ + 1 1 



and their plots are given in Fig. [71 Note that Sj takes a constant value at each complement of 
amoeba, as expected. 

4.3 Topological Charges 



We now move on to discussion of the topological charges given in (12.41) and (12. 5p . For a given 
Laurent polynomial P, these topological charges are evaluated as follows. First let us consider 
the vortex charge. Since the topological charges are independent of the gauge coupling constant 
g, we can take the strong gauge coupling limit g 00. In the strong gauge coupling limit, the 
magnetic flux of the overall U{1) can be written as 



11 1 , , 
Tr F = —ddc log det VL —dd^ log P 

where dc = —i{d — d). By using the Poincare-Lelong formul 



—ddc log \P\ A a 



a, 



X 



X={iu,,U2)e{C*)^ I P(Mi,M2)=0} 



(4.15) 



(4.16) 
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This formula is the generaUzation of the formula dd'^ log \z\ — 2TTS^{z)dx A dy with J dxdy6'^{z) — 1. 
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we can show that the vortex charge can be evaluated as 

V = -c TtFAlu = 2nc tu = 27rcArea(X). (4.17) 

J(C*)2 Jx 

From this computation, it is clear that the vortex charge is distributed on the surface of the 
vortex sheets X. We also see that the vortex charge is uniformly distributed along all X, and the 
total vortex charge is given by the area of the vortex sheets multiplied by the tension 27rc. It is 
interesting to note that the same formula has appeared in mathematics literature ([58], Theorem 
6 and p]). 



We can also give another expression for the vortex charge, using the Ronkin function Np. In 
the tropical limit -Ri, -R2 ~^ 0, the vortex charge is given by 

V = c [ (fx(diTT±i + d2TT±2) = c I d^x{dl + dl) Np{xi,X2), (4.18) 

where c = 4it^RiR2C. Namely, vortex charge is given by an integration of a Laplacian of the 
Ronkin function. If we take the limit Ri = R2 = R ^ 0, the amoeba becomes the tropical 
variety which can be interpreted as the web diagram of the domain walls. The tension of each 
wall can be calculated as follows. Since the integrand of Eq. fl4.18p becomes {df + (9|) Fp{xi, X2) 
in the small radius limit R —>■ 0, the tension of the domain wall is computed by integrating the 
Laplacian of the piece-wise linear function Fp{xi, X2) along the line perpendicular to the wall. If 
the wall is located along the line riiXi + n2X2 + r = {rii +p)xi + (77-2 + q)x2 + r', then the tension 
is given by 

Tip,,) = ^Vp' + q'- (4.19) 

Next let us consider the intersection charge (12. 8p . By taking the strong gauge coupling 
limit, the intersection charge density Xintersection becomes a complex Monge- Ampere measure 
(ddc log |P|)^ on (C*)^ associated with a plurisubharmonicQ function log |P|, which is a higher 
dimensional generalization of the Laplace operator (see [601 161] for discussion on the complex 
Monge- Ampere measure): 

1 1 /9Mog|P|\ 

^intersection = ^^r F A Tr F \ du du- J ^ ^ ^ ^^"^ 

^ ddJog\P\Addc\og\P\. (4.20) 



Then the intersection charge is evaluated again by using Poincare-Lelong formula, 

^ j ddJog\P\AddJog\P\ = dd,\og\Pl (4.21) 



'intersection 



"^^ Monge-Ampere measure is defined for arbitrary plurisubharmonic function. Here it suffices to fcnow that 
fog |P| is plurisubharmonic for arbitrary holomorphic function P. 
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but this naive evaluation is unfortunately divergent. The divergence comes from the fact that 
the strong gauge coupling limit in the master equation (12.141) is ill-defined when f2 = 0, since 
there appears in the master equation. In principle, if we can solve master equation for finite 
gauge coupling, we could safely obtain a correct value of the intersection charge, but that would 
be difficult in practice. Instead, we propose to regularize the divergence as follows. 

Let Pi and P2 be distinguished Laurent polynomials associated with the same Newton poly- 
tope A(P) of P, and replace two P's in (14.211) by Pi and P2, respectively. For generic Laurent 
polynomials Pi and P2, the intersection points of the zero sets of Pi and P2 are discrete points. 
Then we obtain 

^ ^rf41og|Pi| Ad41og|P2| = / d41og|P2| = i#(Xi-X2), (4.22) 



where the surfaces Xi {i = 1,2) are defined by Pi{zi,Z2) = and the number of intersection 
points are denoted as #(Xi ■ X2). 

Thanks to Bernstein's theorem ^3j0 H^{X\ ■ X2) is independent of the choice of the Laurent 
polynomials Pi, P2 as long as Pi and P2 are generic, and is given by 2Area(A). We thus find that 
the intersection charge /intersection is evaluated to be equal to the area of the Newton polygon: 

/intersection = Area(A). (4.23) 

The meaning of this regularization is now clear. The original expression (14.211) is divergent 
essentially because it is a self-intersection number. We propose to regularize this by infinitesimally 
changing P, but with fixed boundary conditions at infinity!^ 

Instead of invoking Bernstein's theorem, we can take more down-to-earth approach and the 
calculation goes as follows. This derivation is not independent from the previous argument and 
moreover not rigorous, but it has an advantage of clarifying the relation with the Ronkin function 
and real Monge- Ampere measure. 

First, it is reasonable to expeclQ that the intersection number does not change under replace- 
ments Pi(mi,U2) Pi(|ui|e*^i, |M2|e*^2) ^nd P2{ui,U2) — »■ P2{\ui\e^'^'^ , |w2|e*'^^), as far as 6'i, 6^2,01 
and 02 are sufficiently generic: 

#(Xi-X2) = #(Xi(0i,^2)-X2(0i,02)), (4.24) 

where 

Xi(^i,^2) = {{z,,Z2)eiC*f\Pii\ui\e''\\u2\e''')=0}, (4.25) 

^2(01, 02) = {{zi,Z2) e (C*)2 I P2{\u^\e''^\\u2\e''f'')=0}. (4.26) 



^"'^This theorem is a generalization of the well-known Bezout's theorem. See [62 for leisurely introduction to 
Bernstein's theorem. 

^^The condition that Pi and P2 are Newton polynomial of the convex polytope A(P) is important. Otherwise 
the answer depends on the choice of Pi and P2- For example, if we take Pi j — zi, P2J = zi + then 
dd'^ log \ Pi, j\ Add'' log \P2j\ = Ofor allj. If we take instead = zi + Z2/j, then dd'^ log \ Pi j\ Add'' log IP^j] = So = 
5'^{zi)S^{z2)dxiAdyiAdx2Ady2- And for ^ Z1 + Z2 (in a neighborhood of 0), dd'^ log |Pij| Add^ log |P4j| ^ jSo- 
All these functions converge to the same P = zi in the j — > cxd limit, but gives a different answer. We thank 
Alexander Rashkovski for providing us with this example. 

^^Essentialy, we are again using Bernstein's theorem here for the rigorous argument. 
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Certainly the intersection number (14.241) might change if 9i and 62 are non-generic, but those 
special values do not contribute when we integrate over all 61 and ^2- The same apphes to 0i 
and 02- Therefore the intersection charge can be written as 

ddc log I -Pi I A ddc log | P2 \ 

,)2 



87r2 
1 [ dOi de2 



I de,M2d^d^ I rf4iog|p,(|^,|e«^M«,|e^^7| Arf41og|P2(ki|e^^M«2|e^^^)| 

1 



Stt 



) 

ddcNp^{xi,X2) A ddcNp^{xi,X2) 



dz, A d^, — — Arp^(xi,a;2) A dz^ A d^/ — — Arp,(a;i, 0:2) 

O/i y J-^2^rp2 f-yXj ^X}^ ^Xl 

d d d d 

dxi A dyi A (ix2 A (i?/2 eikeji — -^Np^{xi,X2) ^ — — A^P2(xi,X2) 



R1R2 [ dx,AdX2flMA{Pl,P2), (4.27) 

7r2 



where 



is known as a real Monge- Ampere measure which is a real analogue of the complex Monge- 
Ampere measure we explained previously. Note that this is well-defined since A'"p is convex, as 
discussed previously. If we set Pi, P2 P, then 

/intersection = -^1-^2 / dxi A dx2 Hessiau (A^p) (4.29) 

= j ci(PiTrSi) A rf(i?2TrS2) (4.30) 
= Area(A(P)). (4.31) 

In the final line we used the fact that PiP2(TrSi, TrS2) takes values in the Newton polytope 
A(P), as explained previously. Interestingly, this result, that the total integral (or "mass" in the 
standard mathematics literature) of the real Monge- Ampere measure is given by Area(A(P)), 
has appeared previously in the mathematics literature ([58] Theorem 4). 

From the above calculation, we have shown that the evaluation of the complex Monge- Ampere 
measure reduces to the evaluation of the real Monge- Ampere measure. Since the Ronkin function 
is described by the zero mode in the KK decomposition, we learn from this fact that contributions 
from KK modes, although present, are canceled out in the final expression. This is just the same 
as in the discussion in Sec. [2) topological charge is determined only from the boundary condition 
and all KK modes vanish at spatial infinity. 

We can provide one more different explanation of the formula /intersection = Area(A(P)), which 
is much easier to understand (although strictly speaking, this is also just a restatement of the 



More precisely, when Pi ^ Pi-, this is called a real mixed Monge- Ampere measure. 
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previous explanation). Taking the limit Ri,R2 oo, the intersection charge /intersection is given 
by the half of the intersection number of the tropical varieties of Pi and P2, and it is easy to see 
that the number is given by 2Area(A(P)) (see Fig. [8] for example). In tropical geometry, this 
statement is known as the tropical Bernstein theorem ([64] Theorem 9.5). 




Fig. 8: Intersection of one tropical variety and its shift in generic directions. The corresponding 
Newton polytope A(P) is given below. It is easy to see that the number of intersection points is 
given by 2Area(A(P)). 



Although the computation here applies only to i?2 ^ cxd limit, we expect that the inter- 
section charge is still given by the same formula for finite Ri,R2 as well, since the intersection 
charge is quantized and does not change continuously depending on Ri,R2. For the same reason, 
although all the arguments so far are in the strong gauge coupling limit g (yo, we expect the 
same formula is kept in the finite gauge coupling as well. 

So far, we have concentrated on the case of Nq = Np = N called local theory. We expect 
that a similar formula holds even for Np > Nq, which is often called semi-local theory. For 
Nq = 1, A^F = 2, we can actually show this rigorously!^ In the strong gauge coupling limit the 
solution of the master equation Eq. (12.141) for Hq = (Pi, P2) is given by 

n ^ ilo = \Pf = \Pi\^ + \P2f . (4.32) 

In this case, /intersection nceds no regularization. We use two representations of the current (inter- 
section charge density) [P] for P = (Pi,P2). One is King's formula (a.k.a. the vector Poincare, 
or the Poincare-Martinelly formula [HIIES]) 

[P] = rfd'=log|P| Arfrf'log|P|, (4.33) 
^^We thank Alexander Rashkovskii for providing us with the following proof. 
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and another is 



[P] = dd' log I Pi I A dd' log IP2I, (4.34) 
which follows from Poincare-Lelong formula applied to P2 on Xi. We thus have 

/intersection = ^ f dd'^Og \P\ A dd'^log \P\ = ^ [ dd'log |Pi| A dd'log jPsj (4.35) 

= Area(A(Pi),A(P2)), (4.36) 

where in the last line we have again used Bernstein's theorem, and Area(A(Pi), A(P2)) is a mixed 
volume defined by 

2Area(A(Pi), A(P2)) = Area(A(Pi) + A(P2)) - Area(A(Pi)) - Area(A(P2)) (4.37) 
and the sum (Minkowski sum) A (Pi) + A(P2) is defined by 

A(Pi) + A(P2) = {x + y\x e A(Pi), y e AiP^)} (4.38) 

For general A'^c and A^f (> A'c), the solution of the master equation Eq. fl2.14p in the strong 
coupling limit is given by 

detn detfio = |Pp = ^|P^|^ (4.39) 

i 

P, = er.,r,-r.^Hl-^^Hl-^- ■ ■ ■ i7o'"«^"« (^ = 1, " " " , A^p l/A^c ! ( A^p - A^c)!). (4.40) 

Therefore the intersection charge /intersection can be evaluated by integrating {ddclog\P\y with 
P being an Np\/Nc\{Np — A"c)!-dimensional vector. Unfortunately, we have no mathematical 
estimate for A^p > 3, and only the upper bound [66] for the integral of the complex Monge- 
Ampere measure in (not (C*)^ is known. We conjecture that the result holds in this more 
general case in the same mannero 

This concludes our discussion of the topological charges. Our discussion mainly concentrates 
on the case A^ = 1. As we have seen, however, even for A^ > 1 case, the overall U{1) part 
represented by Ci = — ^^TrP is still described by the language of the amoeba and tropical 



At least, we can prove that intersection charge /intersection is quantized. We thank Alexander Rashkovskii 
for providing us with this argument. The argument goes as follows. Let P be a iV-dimensional vector P = 
(Pi, P2, . . . , Pjv). For a two-dimensional subspace A of the A^-dimensional complex space (C*)^, let P^ = 
(Pi^,P2^) be P/ = >'ioPi^ where Ay is the 2 x iV-matrix of orthonormal basis of P. Then by [611 [67], 

Monge- Ampere measure is represented as 



{ddMP\f= I {dd,\0g\PA\fd^l{A), 
Jg(N.2) 



IG{N,2 

where dii{A) is the Haar measure on the complex Grassmannian G{N,2). Note (ddc log | Pa |)^ is well-defined for 
all A except for an algebraic subset of G{N, 2) which is zero measure. Since Pa has exactly two components, it 
follows from the discussion of iV = 2 case that J {ddc\og\PA\y is given by the intersection number of P^ and 
P^, which is an integer. This means A 1-^ J {ddc log |Pa|)^ is an integer-valued continuous function. We have now 
proved that the value of /intersection is quantized. 
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geometry. At the same time, we should not think that is the whole story. Out of all the 
topological charges (12.51) . (12.81) and (12. 9p . the instanton number 

/instanton = / = ^ ^ [Tr F A Tr F - Tr (F A F)] (4.41) 

vanishes in U{1) theory and appears only in non-Abelian theory. We will discuss this in the next 
section, but before that let us discuss metric of moduli space for the web of vortices. 



4.4 Metric on Moduli Space 

In this section, we discuss the metric of the moduli space. The metric of the moduh space is 
given by a formula similar to the one-dimensional case, Eq. (13.111) and Eq. (13.121) . 

Kfj = cJ (^^didjlogdetn + 21'^uj AiTi d{ndn-^)dj{ndin-^) - d{ndin~^)dj{ndn-^) ^(4.42) 

This is a new result which is valid for arbitrary values of Nq and Np. 

Here we focus on the case Nq = Np = 1 in which all the moduli parameters are contained 
in the Laurent polynomial. Let us consider a Newton polytope A(P) associated with a Laurent 
polynomial P. Let Vex{Q) and Vin{Q) be sets of external and internal vertices of A(P), respec- 
tively. The coefficients of the Laurent polynomial P are identified with the moduli parameters 
and there exist zero modes corresponding to these moduli parameters. 

PiuuU2)= Yl an„n,u1'u^\ (4.43) 

The coefficients ani,n2! (^i? ^2) G Vex{Q) determine the positions of the external legs and the size 
of the loops of the vortex web, while the coefficients ani,n2! (^i? ""-2) G Vin{Q) determine only the 
size of the loops of the vortex web. Since the zero modes associated to the motion of the external 
legs are non-normalizable, we cannot define the metric for these zero modes and we must fix the 
moduli parameters ani,n2) {'ni,n2) G Vex{Q)- 

First let us consider the case where the loop sizes and the radii of the torus are much larger 
than the width of the vortex sheets / = Xj g^. In this case we can evaluate the leading terms in 
the Kahler metric by taking the thin vortex sheet limit / 

lim Kq = 2c [ (fx didj \og\P\. (4.44) 



1^0 



From Eq. (l4.9p we obtain the Kahler potential in the thin vortex sheet limit / ^ 00 as 

K^87i^cRiR2 J cPx (^Np{xi, X2, a,d) — f{xi, X2, a) — f{xi, X2, a)J . (4.45) 

Here f{xi, X2, a) is a holomorphic function with respect to a„i,n2 5 (^i? ^2) ^ Vin{Q) and it should 
be chosen to make the Kahler potential finite. A possible choice of the function f{xi,X2,a) is 

f{xi,X2) = ^Np{xi,X2), (4.46) 

P(M1,M2) = Yl an,,n2U?uT- (4-47) 

{ni,n2)eVex{Q) 
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Therefore the Kahler potential can be evaluated by integrating the Ronkin functions. 

Next let us consider the case where the loop sizes are much larger than the radius of torus 
R = Ri = R2 ^ I. In this case the leading terms in the Kahler potential can be evaluated by 
taking the small radius limit of Eq. (14 .45 p . Since K ^ 0{R) in the limit R ^ 0, the leading term 
in the Kahler potential is given by 



K ^ STT^cRMm / cfxR (^Np{xi, X2, a, a) — Np{xi, X2)) 

R^O J 



where Fp and Fp are piece-wise linear functions defined by 



Fp 



lim RNp(xi,X2) 

lim RNp(xi,X2) 
R^o ^ 



max [riiXi + ^2X2 + Tm.no), 
max {niXi + n2X2 + Tn^^nJ, 

{n\,n2)&Vex{Q) 



(4.48) 

(4.49) 
(4.50) 



with Tr, 



,rt2 — ^ I ,n2 \ ■ 

Let us consider the simplest example of one-loop configuration associated with the Laurent 
polynomial 



P{Ui, U2) = Ui + U2 + ^U2 ^ 

In this case the vertices of the Newton polytope are 



(4.51) 





ao,i) 


(-1,-1) 


(0,0^ 


^1,0) 




(a) Newton polytope (b) web diagram (tropical variety) 



Fig. 9: Newton polytope and web diagram for P = ui + U2 + ^ + ao,o- The loop consists of 

three vortices, which can be interpreted as walls with tension Ti = T2 = c/ R and T3 = \/2c/R, 
see Eq. flCTj) 



V{Q) = {(-1,-1), (0,0), (1,0), (0,1)} (4.52) 
V^niQ) = {(0,0)} (4.53) 
Ve.{Q) = {(-1,-1), (1,0), (0,1)} (4.54) 



24 



There exists only one normalizable moduli parameter ao,o for which we can define the metric. 
This moduli parameter ao,o is related to the size of the loop in Fig.[9]-(b), which is proportional 
to r = i?log |ao,o|- The non-normalizable moduli a_i^_i, ai^Q, ao,i have already been fixed to 1 in 
Eq. ( I4.5ip . The Laurent polynomial defined in Eq. (I4.47P is given by 

P{Ui,U2) = Ui + U2 + U^^U2^. (4.55) 

The piece-wise hnear functions Fp{xi,X2) and Fp{xi,X2) defined in Eq. ( ]4.50p are now given by 

Fp(xi, X2) = max(-xi - X2, xi, X2, r), Fp{xi,X2) = max(-a;i - X2, xi, X2), (4.56) 

By integrating Eq. (14.481) . we obtain the asymptotic Kahler potential as the volume of a tetrahe- 
dron in Fig.[10l which is given by 

K^Un^cRr^ (4.57) 
Differentiating the Kahler potential Eq. ( 14.57^ . we obtain the metric of the moduli space and 




Fig. 10: The asymptotic Kahler potential is proportional to the volume of the tetrahedron 
surrounded by four planes f{xi,X2) = —Xi — X2, f{xi,X2) = Xi, f{xi,X2) = X2, f{xi,X2) = r. 

then the effective Lagrangian 

Leff = ISTT^crR (^r^ + R^O^^ . (4.58) 

Here 9 is the phase of the moduli parameter a = e*"/^"^*^. This asymptotic Lagrangian can be 
interpreted as the kinetic energy associated with the motion of the three vortices composing the 
loop. Since they can be interpreted as walls with tension Ti = T2 = c/R and T3 = \/2c/ R, 
two of three vortices have the mass mi^2 = Ti^2 x length = 3cr/R and the other vortex has the 
mass 7713 = ^3 X length = 6cr/R. If the moduli parameter varies with velocity r, the vortices 
move in the (xi, X2)-plane with velocities Vi^2 = 'i^ and ^3 = r/^/2. Then the total kinetic energy 
associated with the motion of the vortices is Yl^=i ^'^1 — ^cr/Rr"^ = ISir'^crRr'^. 

5 Instantons inside Non-Abelian Vortex Webs 

In this section, we discuss the instanton number /instanton = / C2 as promised in the previous 
section. We consider the case of U{2) gauge theory with Np = 2 scalar fields as the simplest 
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model which admits the BPS configuration with /instanton 7^ 0. The generahzation to the U{N) 
gauge group should be straightforward. 



5.1 Instanton Number on a Planar Vortex Plane: a Review 

Let us first review the 1/2 BPS vortex moduli space. The vacuum of U{2) gauge theory 
with A^F = 2 scalar fields breaks the color and fiavor symmetry U{2)c x SU(2)p into color- 
fiavor locked symmetry SU{2)c+f- The 1/2 BPS single vortex in this theory further breaks 
SU{2)c+F into f/(l)c+F- Therefore there appear Nambu-Goldstone modes of the complex pro- 
jective space CP^ ~ SU (2) c+f/U{1)c+f localized around the vortex [21 [3]. The orientation 
moduli orientation = CP^ form a part of the moduli space. 

The moduli matrix for the vortex at 2:1 = with an orientational moduli parameter b is given 
by [5] 

where ~ represents the V^-equivalence relation. These two moduli matrices provide two patches 
b and 1/6 of CP^ [5]. For this moduli matrix, the solution of the master equation is given by 

where is a solution of the master equation in the case of the Abelian-Higgs model {Nq = 
Np = 1), such that 

d,,{n,d,,n;') = -^(i - \z,\'n:'). (5.3) 

Therefore once the vortex solution Q^, of the Abelian-Higgs model is given, one can construct the 
whole solution of the non-Abelian model. 

Next let us discuss the 1/4 BPS configuration of instantons and vortices on (C*)^. The 1/2 
BPS vortex plays the role of a host soliton in the 1/4 BPS vortex-instanton configuration. It 
has an internal degree of freedom parametrized by an orientational moduli parameter b, which 
is interpreted as an inhomogeneous coordinate of A^orientation = CP^ as denoted above. The 
instantons inside the vortex can be constructed as lumps in the vortex effective theory which 
is the CP^ sigma model [221 123]. By using the moduli matrix (15.11) for 1/2 BPS single vortex, 
we can construct the moduli matrix for some 1/4 BPS configurations as follows. The moduli 
matrix for the 1/4 BPS vortex-instanton configuration is obtained by promoting the orientational 
moduli 6 to a holomorphic function of the other holomorphic coordinate Z2 [23] 

i^o = v^(j ^^'^fY b{z,) = af[{z, - z!i\ (5.4) 
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where the moduh parameter ^g*"* denotes the position of the i-th instanton and a determines 
the overall size of the instantons on the vortex sheet defined by Zi = 0. This function b{z2) is 
regarded as a holomorphic map from the vortex sheet to orientation = CP^. 

In order to compute the instanton number, we require the information upon the solution 
of the master equation Q for the moduli matrix Eq. fl5.4p . Since the topological charges are 
determined only from the boundary condition, it can be calculated from the asymptotic behavior 
of the solution Q which is obtained by the following procedure. Solving the BPS equation in 
the 2+1-dimensional effective theory on the vortex worldvolume, we obtain the lump solution 
6(22) corresponding to the instanton inside the vortex sheet. Substituting back the solution of 
the effective theory b{z2) into the 1/2 BPS vortex solution Eq. (15. 2p . we find 

Although this is not an exact solution of the master equation since h{z2) is not a constant, 
we can see that this matrix Vl possesses the correct asymptotic behavior at spatial infinity by 
substituting Vl into the master equation. Therefore the topological charge, which is determined 
by behavior of the fields at spatial infinity, can be evaluated from this matrix VL. Inserting 
F = -iS^^d{ytdn-^)S into Eq. (El, we obtain 

^ ^Tr(FAF) = -i^ /d41oga(2i,^i)Arf41og(l + |6(z2)r) = fc. (5.6) 



^7r2 / ' ' 87r2 



Therefore the instanton charge is measured by the degree of the function 6(2:2); namely the degree 
of the holomorphic map from the vortex sheet C* to A^orientation = CP^. 



5.2 Instanton Number on Non-Abelian Vortex Webs 

As we have seen in the example above, the instantons can exist on the planar vortex sheets. Now 
we work out more general 1/4 BPS configurations explicitly. The planar host vortex sheet can 
be extended to a general web of the vortex sheets characterized by P{ui, U2) as discussed in the 
previous sections. We therefore replace the lower-right component zi in the moduli matrix (15.41) 
by P{ui,U2). Moreover the holomorphic function b in (15. 4p is also replaced with a function of 
two holomorphic coordinates Zi and Z2- The moduli matrix for such 1/4 BPS configuration of 
the instantons and vortex sheets on (C*)^ thus becomes 

V P{Ui,U2) J 

where P{ui,U2) and b{ui,U2) are Laurent polynomials 

P{Ui,U2) = ^ttn^^niUTuT, b{Ui,U2) = bn^ ^uiUT ■ (5-8) 

Although Eq. (15. 7p is not the most general form of the moduli matrix for the 1/4 BPS configura- 
tions, we treat this simple form of the moduli matrix for essential explanation. The coefficients 
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of the Laurent polynomials a^^^nj ^ni,n2 are the moduli parameters which give the location 
of vortex sheets, and the positions and sizes of the instantons, respectively. 

For the moduli matrix Eq. (15.71) . the instanton charge / in Eq. (12. 4p is computed as follows. 
Since the topological charge should not change under continuous deformations, we can take a 
strong coupling limit (7 — oo in which the solution of the master equation Q approaches to 

. ^ ^ ^H^Hl ^ ( 1 ;f |5 ) . (5^9) 

However a direct calculation does not work since the vortex sheets become singular in this limit. 
To avoid the calculation involving the singular vortex sheets, we perform another deformation of 
the configuration. Similarly to the case of the planar vortex sheet, let be a 2 x 2 matrix given 
by 

where is a solution to the following equation 

2 

d,,{n,d,,n;')+d,,{n,d,,n;') = -^(i - \p\'n:'). (5.11) 

Although the matrix Q is not a solution of the master equation unless the holomorphic function 
b{ui,U2) is constant everywhere, this matrix has the correct topological information about the 
configuration similarly to the planar case. For this matrix Q, we can show that 



J ch2 = j Ti [d{ndn-^) A d{ndn-^)] 



„ , . ddJogQ^ A ddJog{l + \bf) - l-ddJogQ^ A ddJognA . (5.12) 

lOTT^ J \ 2 J 

If we take the strong gauge coupling limit — >■ oo in Eq. (15.111) . then approaches to |-Pp. 
Therefore we obtain the instanton charge / in Eq. (12. 4p as 

^=^j {ddclog\P\Addclog{l + \b\'')-ddclog\P\Addclog\P\), (5.13) 

where the first term gives the instanton number /instanton in Eq- (12.91) and the second term gives 
the intersection charge /intersection in Eq. (12.81) which has been computed in Sec. 14. 31 The instanton 
number is rewritten by using the Poincare-Lelong formula 



'instanton 



^[ dd Jog\P\ A dd Jog{l + \b\') = ^ [ ddJog{l + \b\^), (5.14) 
7^ J{C')2 47r Jx 



where X denotes the zero locus of P corresponding to the vortex sheets. Therefore the instanton 
number is given by the degree of the map b\x : X ^ CP^. 

To see a distribution of the topological charge density, we take two limits of the parameters: 
one is a small instanton limit and the other is a small radius limit /? — > 0. The small instanton 
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limit is realized by taking the limit &ni,n2 ~^ ^ with fixed ratios &ni,n2/^ni,ri2- Then the two form 
ddclog{l + ddc\og\b\'^ has a delta function-like support on the zeros of h{ui,U2)- From 

this fact we find that in the small instanton limit the instantons are localized at common zeros 
of b{ui,U2) and P{ui,U2), and that the vortex sheets are located at P{ui,U2) = 0. We also find 
that instantons are localized at the positions of lumps from the viewpoint of effective theory on 
the vortex sheets. 






I = 1 
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/ = 1 











(b) 











/ - 1 









U = 2 



(d) 



(e) 



Fig. 11: The instanton number density in the small radius limit for P = Ui + U2 + 1 and 
b = 6i,iMiM2 + bifiUi + 6o,o- The instanton number density is localized at the intersections of the 
tropical variety for P (solid line) and the lines on which Fh{xi,X2) is not differentiable (dashed 
lines) . 



Next, let us consider the small radius limit i? — > 0. In this limit, the function log(l + |6p) 
becomes 

R 

— log(l + Fh{xi, X2) = max (niXi + 112X2 + Sn.^^J, (5.15) 

Z (ni,n2) 

where So,o = "f log(l + l^o.oP) and Sni,n2 = -Rlog I^ni,n2l ^01 (ni,n2) 7^ are fixed in the limit. 
Then the instanton number takes the form 



' instanton 



87r2 



d4iog|p| Ad4iog(i + i&r) 

(C*)2 

d^xeijeki^^—Fp{xi,X2)^—^Fb{xi,X2). (5.16) 
OXi OXk OXj 0x1 

Therefore the instanton number density is localized at the intersection of the tropical variety of 
the Laurent polynomial P and the lines on which the piece- wise linear function Fb{xi,X2) is not 
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different iable. Wlien tlie instanton number density is localized at the intersection of the lines 
riiXi + n2X2 + r = and niXi + n2X2 + s = 0, the instanton number is given by 

/instanton = j (fxeijeki ninkS(niXi + n2X2 + r) hjhi5{niXi + ^2X2 + s) 

= \n1n2 - nin2 \ . (5-17) 



det 



Til n2 
fii h2 



Fig.[TT] shows examples of the instanton number density on for P = ui + U2 + 1 and h = 
6i_iMiM2 + 61,0^1 +&o,o- Varying the moduli parameters, the instantons localized at the intersection 
move along the tropical variety of P. For each intersection in Fig. [T]T a)-(d) the instanton number 
is / = 1. In Fig. fTTT e) 1 = 2 instanton, which can be interpreted as the coincident instantons, is 
localized at the intersection of the lines xi — X2 = and xi + X2 + -Rlog — ^ log(l + |6o,oP) = 0. 

So far we have treated the specific configuration of the non-Abelian vortex sheets and instan- 
tons on (C*)^ and the specific gauge group U{2). For more general cases, the computation of the 
instanton number seems to be complicated and difficult, but the essence should be similar to the 
above calculations. 



6 Conclusion and Discussion 

In this paper, we have investigated generic intersections (or webs) of vortices with instantons 
inside, which is a 1/4 BPS state in the Higgs phase of (the bosonic part of ) five- dimensional 
Af = 1 supersymmetric U{Nc) gauge theory on Rt x (C*)^ ~ M^'^ x with A^f = Nc Higgs 
scalars in the fundamental representation. We have found that this vortex-instanton system 
can be beautifully and naturally understood in the mathematical framework of the amoeba and 
tropical geometry, and have proposed a dictionary relating the solitons and gauge theory to the 
amoeba and tropical geometry (summarized in Tabled]). 



Table 1: Dictionary relating soliton/gauge theory to amoeba/tropical geometry. 



soliton/gauge theory 


amoeba/tropical geometry 


moduli matrix Ho{zi, Z2) 
projection of vortex sheet 

position of step-wise kinks 
Wilson loop TrSj 
intersection charge /intersection 
vortex charge density V 


Newton Polynomial P{ui,U2) 
amoeba Ap 
tropical limit 
tropical variety 
derivative of Ronkin function: diNp, f|4.9|) 
(total mass of) complex Monge- Ampere measure 
Laplacian of Ronkin function: fl4.18p 



In this discussion, the moduli matrix formalism has played crucial roles. The solutions to 
1/4 BPS equations are parametrized by a holomorphic function (Laurent polynomial) Hq{zi,Z2) 
of two complex parameters zi, Z2 of (C*)^. This Laurent polynomial can also be considered as a 
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Newton polynomial of some convex polytope A, namely the grid diagram. In the strong gauge 
coupling limit, the position of vortices is exactly given by the zero of Hq^ while the projection of 
the shape of vortex sheet is the amoeba of A. Moreover, we can relate Wilson loops in T^, or 
the zero modes of gauge fields in Kaluza-Klein decomposition, to the derivatives of the Ronkin 
function Nhq{xi, X2), which is a convex function and is defined from Newton polynomial Hq. 

The relation with the tropical limit and tropical geometry has also been discussed. In the 
discussion of solitons, it is natural to consider dimensional reduction of the theory in order to 
obtain BPS solitons in lower- dimensional field theories. This limit is known in the mathematical 
literature as the tropical limit. In this limit, the shape of amoeba degenerates into a tropical 
variety, which is nothing but the so-called {p, g)-web of the grid diagram. We have shown that 
tropical geometry provides simple and elegant method to understand not only the dimensionally 
reduced theory but also the original vortex-instanton system on (C*)^. 

We have also discussed the topological charges, which are divided into three types. They are 
the vortex charge V, the intersection charge /intersection (negative contribution of the instanton 
charge) and the instanton number /instanton- First, the vortex charge is uniformly distributed 
along vortex sheet det//o = 0, and its density is given by the Laplacian of the Ronkin function in 
the tropical limit. Its total charge is the area of the vortex sheet multiplied by 27cc. Second, in 
the strong gauge coupling limit the intersection charge density is given by the complex Monge- 
Ampere measure of a plurisubharmonic function log \Ho{zi, Z2)\, and total intersection charge is 
given by the area of the grid diagram A with a suitable regular izat ion. Third, the instanton 
number /instanton appears only in the non-Abelian case, and we have discussed the case of Np = 
Nc = 2 as an example. Our discussion simplifies in two limits. In the small instanton limit, 
instantons are localized at intersections of P{zi,Z2) and another Laurent polynomial b{zi,Z2), 
which parametrizes an orientational moduli CP^. In the small radius limit, the instanton number 
density Xinstanton is localized at the intersection of the tropical varieties corresponding to P and 
b. We have also obtained the general form of the Kahler potential and the asymptotic metric of 
the moduli space of a vortex loop as a byproduct of the discussions above. In the tropical limit, 
the Kahler potential is given by the volume of a convex polytope, and the effective Lagrangian 
can be interpreted as the kinetic energy associated with the motion of vortices composing the 
loop. 

In our discussion, we mainly focused on the case of Abelian-Higgs model {Np = Nq = 1), 
but as far as the overall U{1) part is concerned the story is exactly the same in non-Abelian 
U{Nc) case. We also obtained new results by going to non-Abehan gauge group, such as the 
instanton number. This seems to suggest a non-Abelian generalization of the amoeba and tropical 
geometry. 

Despite such impressive success, there are still many points which need further exploration. 

First, in this paper, we have found one-to-one correspondence between the amoeba/tropical 
geometry and solitons in Abelian gauge theory. However the configuration of instantons inside 
non-Abelian vortex- webs in non-Abelian gauge theory discussed in Sec. Odoes not correspond to 
the amoeba and tropical geometry so far. This configuration suggests non-Abelian generalization 
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of amoeba and tropical geometry. Furthermore the non-Abehan vortices have been recently 
extended to the case of gauge group G = U{1) x G' with G' arbitrary simple group [12]. This is 
the case of complex one dimension. It should be extendible to the case of complex two dimensions 
such as (C*)^. That may suggest further generalization of the amoeba and tropical geometry 
associated to arbitrary group, which contains the usual one as a special case of U{1) gauge 
group Ezl 

Second, let us pursue a possibility to generalize the space where the amoeba lives. In this paper 
we have considered the amoeba on (C*)^. This is because the maximal space-time dimension of 
supersymmetric gauge theory with eight supercharges is (i = 5 + 1. Since space-dimension 5 is 
odd we have studied four dimensional case (C*)^. If we abandon implementing supersymmetry, 
we can extend the bosonic Lagrangian (12. ip to space dimensions higher than 5, and study higher 
co-dimensional composite solitons of the vortices and instantons extending to various directions. 
In fact the generalized vortex equations on arbitrary Kahler manifold of arbitrary dimensions 
were obtained in (bosonic) Yang-Mills-Higgs theory |39j- Therefore we expect those equations 
give further correspondence of the amoeba and gauge theory on (C*)" or on general Kahler 
manifolds. 

The third topic is the relation with dimer model [68l [69] . Dimers do not appear directly in our 
discussion, but it is known that dimer model is intimately connected with the amoeba and tropical 
geometry. For example, the Ronkin function as defined in (14.91) coincides with the thermodynamic 
limit of partition function of a dimer model. Moreover, a spectral curve of the dimer model is 
known [68] to parametrize Harnack curve, whose amoeba has good properties. The dimers also 
appear in discussion of the brane tihngs [70], [7^ and four- dimensional J\f = 1 superconformal 
quiver gauge theories. The brane tilings now have an interpretation as configuration of D5- 
branes and NS5-brane [72], whose brane configuration is shown in Table [21 wherein S is a two- 
Table 2: The five-brane configuration described by brane tilings. In the weak string coupling 
limit, the surface S becomes a zero locus of a Newton polynomial corresponding to the toric 
diagram. This setup is analogous to our vortex-instanton systems, although details are different. 








1 2 


3 


4 5 6 7 


8 9 


D5 


o 


o o 


o 


o o 




NS5 


o 


o o 


o 


S (2-dim surface) 





dimensional surface in 4567-directions. In the weak gauge coupling limit, this S is the zero locus 
of a Newton polynomial with respect to two complex variables in (C*)^. Hence the mathematical 
structure of the NS5-brane is exactly the same as that of the vortex sheets we have considered in 
this paper. The toric diagram and the grid diagram are identical, and the meaning of the tropical 
limit and tropical variety also coincides. Of course, we should keep in mind that there exists 

^^In fact, the motto of tropical geometry is to extend usual algebraic geometry by replacing commutative ring 
with a commutative semiring. Another generalization is to replace a commutative ring with a non-commutative 
ring, which is non-commutative geometry. Perhaps our discussion of non-Abelian vortices suggests further gener- 
alization by combining above two, which should be called "non-commutative/non-Abelian tropical geometry". 
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important differences between the soliton systems in this paper and the brane tihngs. First, 
in this paper we have assumed systems with eight supercharges, but in brane tihng, we have 
only four supercharges. Related to this fact is that we do not have analogue of D5-brane in the 
soliton side, and have the instanton charge instead. Still, we might obtain something new from 
this analogy. For example, in the discussion of brane tilings, the projection of (C*)^ on to 
directions, which is called alga or coamoeba, plays crucial roles [371 ES]. It would be interesting to 
see whether coamoeba has any significance in our setup. Perhaps we can understand these points 
better if we can find a D-brane realization of the vortex-instanton system. For this direction, 
the work on D-brane configuration of vortices on cylinder [71 |TT], its T-dual to the D-brane 
configuration of the domain walls [151 the D-brane configurations of the domain wall 

webs [28] should be a useful guideline. 

The amoeba and tropical geometry also appear in the computation of the topological string 
amplitude and the instanton counting in 4d and 5d supersymmetric gauge theories. Their per- 
turbative dynamics are ruled by asymptotic behavior of the (plane) partitions where the amoeba 
and the Ronkin function appear. It is also pointed out that the Kahler structure behind the 
theories is closely related with the volume of the convex cone of the Ronkin function similarly 
to our discussion on the Kahler potential. These relations suggest that the intersecting soliton 
system also admits an interpretation of the microscopic partitions (dimers) and that the Kahler 
geometry of the moduli space of the solitons is determined by a suitable asymptotic limit of the 
microscopic interpretation. 

Finally, statistical partition functions of the vortices on a cylinder were studied by using D- 
brane configurations [TT]. There the integration over the moduli space of vortices is drastically 
simplified in the T-dual picture: vortices are mapped to domain walls and the integration reduces 
to a problem of rods. The limit of parameters employed in [TTj is a little different from that in 
this paper. There the limit g ^ oo, R ^ was taken with fixed wall width d = while 
in this paper we have taken the strong coupling limit g ^ oo first and then the small radius 
limit i? — 0. This is nothing but the non-linear sigma model version of the limit in [TT]. In 
the case of Np > Nq, the gauged linear sigma model reduces to the non- linear sigma model and 
the vortex solution reduces to the lump solution [10]. The lump solution on the cylinder can be 
mapped to kinks on one-dimensional space. The small radius limit i? — > 0, which corresponds to 
dequantization (ultradiscretization) limit limji^Q Rlog{e'^^^ + e^^^) = max{A,B), enable us to 
identify the kinks with free particles in one- dimensional space as in Sec. [31 The partition function 
of the free particles gives the exact volume of the moduli space of the sigma-model lumps, namely 
the partition function for a multi-lump system. This fact and the result of [H] suggest that the 
procedure of the dequantization is powerful enough to give the exact partition function of the 
solitons in non-linear sigma model and variant of it makes the computation of the partition 
function very simple even in the case of finite gauge couplings. This should be extendible to 
the vortex webs discussed in this paper, to obtain a partition function of the instanton-vortex 
system. We expect that it will be reduced to the Nekrasov's partition function in the limit of 
g^/c —>■ where the vortices disappear while the instantons still remain. The integration over 
the moduli space of the instanton-vortex system may provide a systematic method to compute 
the symplectic Gromov-Witten invariant, which is a combination of the Donaldson invariant and 
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the Gromov-Witten invariant 
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